Dynamic behavior of anisotropic non-equilibrium driving lattice gases. 
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It is shown that intrinsically anisotropic non-equilibrium systems relaxing by a dynamic process 
exhibit universal critical behavior during their evolution toward non-equilibrium stationary states. 
An anisotropic scaling anzats for the dynamics is proposed and tested numerically. Relevant critical 
exponents can be evaluated self-consistently using both the short- and long-time dynamics frame- 
works. The obtained results allow us to clarify a long-standing controversy about the theoretical 
description, the universality and the origin of the anisotropy of driven diffusive systems, showing 
that the standard field theory does not hold and supporting a recently proposed alternative theory. 
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The development of a theoretical framework for the 
study of non-equilibrium systems is of vital importance 
for a wide range of disciplines in both science and tech- 
nology. In order to overcome such a difficult problem, it 
is useful to study models capable of capturing the essen- 
tial physical features of the real system. The statistical 
mechanics of interacting lattice gases, driven into non- 
equilibrium steady states (NESS) by external fields jiJU, 
has attracted growing attention. Some remarkable char- 
acteristics of the NESS observed in driven diffusive sys- 
tems (DDS) are, among others, their non-Hamiltonian 
nature, the violation of the fluctuation dissipation the- 
orem and the fact that the steady-state distribution is 
determined by the dynamics. 

The prototypical DDS has early been proposed by 
Katz, Lebowitz and Spohn (KLS) |3| as a generaliza- 
tion of the Ising model. Then, the KLS model is a ki- 
netic Ising model with conserved dynamics. Using the 
lattice gas language, particle hopping in the direction 
(against the direction) of an externally applied field (E) 
is favored (unfavored) while jumps perpendicular 

to the field are unaffected by it. For half- filled lattices 
and high enough temperatures the KLS model exhibits 
a disordered phase. However, at low temperatures an 
anisotropic ordered NESS emerges and is characterized 
by strips of high particle density crossing the lattice in 
the direction parallel to the applied drive So, at 

some critical temperature (T c ) the KLS model undergoes 
a second-order non-equilibrium phase transition. The is- 
sue of the universality class of the KLS model is of general 
interest because it will contribute to the rationalization of 
non-equilibrium critical phenomena. However, this task 
is elusive and has become the subject of a long-standing 
debate. A detailed discussion of such conceptually rich 
debate is beyond the aim of this work. Instead, we will 
briefly comment their theoretical and numerical aspects. 

The theoretical controversy . Aimed to analyze the 
critical nature of DDS and determine their universality, 
a Langevin equation has been proposed This equa- 
tion describes the stochastic evolution of the local parti- 
cle density p(x, t) and in terms of 4> = 2p — 1, it reads 



-Vl0 + r ± Vi</»+|A^ 
+T||V|(/»-eV|| ( /» 2 + C(r,i), 



(1) 



where £(r, t) is a conserved noise term that reflects the 
fast degrees of freedom, e denotes the coarse-grained 
drive, and tj_, T|| and g are model parameters. This field 
theory leads to anisotropic scaling wave vectors in the 
critical region, with non-trivial correlation length expo- 
nents vj_ and The current term (eV||</> ) is the most 
relevant nonlinearity and source of anisotropy. 

Another Langevin equation, which is a coarse-grained 
version of the master equation, has also been proposed 
H . Although such equation has been criticized || , after 
healing some flaws j?], the new equation, which is valid 
in the limit of an infinite driving field, reads 



dr 



(r) = -Vl> + r ± Vi0 
+7 1 |Vjto + C(r,t). 



-A i 



(2) 



Comparing Eqs.(l) and (2) it follows that they only dif- 
fer in the current term eV||(/> 2 that vanishes in (Eq.(2)). 
This difference has deep consequences, including: i) The 
critical dimension of Eq.(2) is d = 3 instead of d = 5 
in Eq.(I). ii) Renormalization group analyses give dif- 
ferent critical exponents (see Table I), iii) It has been 
argued that the most relevant non-equilibrium effect is 
the anisotropy generated by the driving field g] and not 
the current Q. iv) In Eq.(2), the current term vanishes 
in the infinite driving limit, as shown rigorously in ||. 

Another piece of heavy debate is the disagreement 
in the interpretation and evaluation of numerical data. 
Valles and Marro |J, have concluded that j3 seemed to 
be close to j3 s=s 0.3, as predicted by Eq.(2). In view 
of the anisotropic nature of the system, Leung ]To| ] has 
developed an anisotropic finite-size approach. Extensive 
computer simulation [10 11 results give data collapsing 



when the shape of the lattice (S = Ly X /L x ) and (3 are 
taken in accordance with Eq.(l). However, a subsequent 
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analysis of Leung's data [JlO| , discussed in , shows that 
they may also be consistent with (3 close to 0.3. Also 
recent simulations give (3 ~ 0.33 S. 

We note that the analysis of numerical data obtained 
under NESS conditions may be biased by some assump- 
tions such as anisotropic vs isotropic scaling, the shape 
S of the lattice, the method used to estimate T c in the 
thermodynamic limit, the quality of data collapsing, etc. 

The aim of this work is to provide an unambiguous 
clarification of the existing controversy, based on exten- 
sive numerical studies of both the short- and the long- 
time dynamics of four variants of the KLS model. In 
their work, Janssen et al. |l^] have shown that a system 
relaxing by a dynamic process exhibits universal behavior 
at the early stages of its evolution towards equilibrium 
states. This idea has been verified in many models ex- 
hibiting critical behavior under equilibrium conditions 
p3| . Since the extension of this concept to the field of 
non-equilibrium critical behavior and to anisotropic 
systems has not been performed yet, the present work 
has additional ingredients of general interest. 

The KLS model j3| is defined on the square lattice 
assuming a rectangular geometry L X ,L V , using periodic 
boundary conditions. A lattice configuration is speci- 
fied by the set of occupation numbers m,j = {0, 1}, cor- 
responding to each site of coordinates (i,j). Nearest- 
neighbor (NN) attraction (J > 0) is considered. So, in 
the absence of a field the Hamiltonian H is given by 



n 



-4J 



n i,j n i',j'i 

<ij,i'f> 



(3) 



where the summation runs over NN sites only. The driv- 
ing field E acts along the y— direction. The coupling 
to a thermal bath at temperature T and the action of 
the field are considered through Metropolis rates, i.e. 
min[l, exp — ({ATL — r/E} /ksT), where ks is the Boltz- 
mann constant, AH. is the change in 7i after the ex- 
change, and rj = (—1,0,1) for a particle attempting to 
hop (against, orthogonally, along) the field, respectively. 

Monte Carlo simulations are performed using lattices 
of different sizes, with 240 < L x < 960 and 30 < L y < 
480. T is reported in units of J/fcs and E is given in 
units of J. In all cases the density is p ~ 1/2. One 
Monte Carlo time step (mcs) involves L x x L y trials. 

Four variants of the KLS model are studied, namely: i) 
the infinite driving limit (IKLS) with E — oo, ii) a finite 
driving case (FKLS) with E = 0.5, Hi) the random field 
KLS model |14|] in the infinite driving limit (IRKLS) and 
iv) the oscillatory KLS model |l5| in the infinite driv- 
ing limit (IOKLS). In the IRKLS (IOKLS), the driven 
field takes values E — ±oo at random (with a period t ), 
generating anisotropy but not an overall current. So, the 
study of the IKLS model is aimed to determine the valid- 
ity of either Eq.(l) or Eq.(2). Also, the FKLS is studied 
to check the validity of Eq.(l) in the finite driving case. 
Finally, the comparison of the results of the IKLS model 
and those of both the IRKLS and the IOKLS models will 
help to clarify the role of the current. 



In order to observe universality in the dynamics the 
system has to be initialized with configurations far from 
the NESS. For this purpose two different initial states 
have been employed, namely: i) Fully disordered con- 
figurations (FDC) as expected for T — > oo and ii) The 
ground state configuration (GSC) as expected at T = 0, 
which is a single strip parallel to the drive |Tq ]. 

The order parameter (OP) is defined as the excess den- 
sity in the direction parallel to the applied field, namely 



OP = (RL X ) 



Po\ 



(4) 



where P(i) = (Ly) -1 Ylf=i n ij ^ s the density profile 
along the x— direction (perpendicular to the drive) and 
R = (2p a (l — p )) is a normalization constant. For sim- 
ulations starting from the GSC the OP given by Eq.(4) 
and that proposed by Leung ||Io|| give the same results. 
However, in contrast to Leung's OP, that given by Eq. 
(4) is suitable to detect the onset of multi-stripped or- 
dering during the short-time dynamics when simulations 
start from the FDC. Since OP involves the calculation 
of the absolute value of the excess density, it is not suit- 
able for the observation of the expected critical initial 
increase in the order parameter fl3|| . Instead, OP can 
be identified with the square root of the second moment 
of the excess density. So, using FDC starting configura- 
tions, the proposed scaling anzats, which generalizes the 
standard one [fl3|, reads: 



OP(t, r, L yi L x ) = b- p/u WP*(b- z t,b 1/u \W, 

b- l Ly,b- V ^L a \ 



(•5) 



where r = T — T c /T c , OP* is a scaling function, b is the 
spatial rescaling factor, (3 is the order parameter critical 
exponent, v\\ (v±) is the correlation length exponent in 
the direction parallel (perpendicular) to the drive, and 
z is the dynamic critical exponent, respectively. In or- 
der to generate the FDC the lattice is filled at random 
with probability p = 1/2. While a procedure is used to 
ensure that the whole density of the sample is exactly 

— 1/2 

p = 1/2, density fluctuations of the order of L y are 
present along the columns parallel to the drive. These 
tiny fluctuations add up according to Eq.(4), so that the 



amplitude of the OP depends on L y as L y 
b cx t 1 ^ in Eq.(5), it follows 

l l u \\ z n 



1/2 



Taking 



OP(t,L y ) cx L y 2 t C2 OP**(t i/v " z T),L x ,L y -> oo , (6) 



with c 2 = (1 - 2j3/v\\)/2z 1 13 1 



Figure 1(a) shows log- log plots of OP(t)L\/ 2 vs t for 
the IKLS model obtained at T = 3.20. It should be 
noticed that a tiny downward (upward) deviation from 
linearity is obtained for T = 3.21 (T = 3.19) (not shown 
here). So, our estimation of the critical temperature is 
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T IKLS = 3 20 ± o Q1 g0; T IKLS _ L41T / 5 where T I 

is the Onsager temperature of the d = 2 Ising model. 
Notice that this value is fully consistent with previous 
estimations J|,[|[l0]|ll[] . Data collapsing of the raw data 
already observed for different lattice sizes supports the 
anzats of Eq.(6), and this result confirms the assump- 
tion, also implicit in Eq.(6), that for the case of FDC the 
only relevant length of the lattice is that parallel to the 
drive along which the precursors of the striped patterns 
start to develop. The best fit of the data shown in Figure 
1(a) gives c 2 = 0.114 ± 0.005. Notice that all the critical 
exponents evaluated in this work are listed in Table I. 



PLh' 

o 



• 30x240 




+ 60x240 




□ 90x240 




120x240 




» 240x240 




♦ 120x30 


/ 


■ 120x60 




▼ 120x120 














(a) 




10 10 10 no 



FIG. 1. Log-log plots of the scaled order parameter vs t 
at criticality. Results averaged over 10 4 different realizations, 
(a) Data corresponding to the IKLS model obtained starting 
with FDC and using lattices of different sizes as shown in the 
figure, (b) Data corresponding to the four studied models 
and obtained using lattices of size L x — 240, L y — 120. The 
straight line has been drawn to show that after some tran- 
sient, the curve of the FKLS model becomes almost parallel 
to those of the other cases. 

Figure 1(b) allows the comparison of data obtained 
for the four different models studied. A perfect over- 
lap of the data corresponding to the IKLS, the IRKLS 
( T irkls = 3 16 ± o 02) and the OKLS (T° KLS = 
3.16 ± 0.02) is found, leading to almost the same expo- 
nent C2, as shown in Table I. Data of the FKLS model 
( T fkls = 2J8±0.015, E = 0.5) exhibit a crossover from 
a very early behavior with appreciable curvature to a lat- 
ter regime (t > 10 3 ) where the slope of the other cases is 
recovered (within error bars, see Table I). 

Furthermore, taking the logarithmic derivative of OP 
given by Eq. (6) at criticality, one has |]l3| 



dlnOP(t, t) 



(xt c 



(7) 



The best fits of the data give the values of the exponents 
C|| listed in Table I. 

Starting from GSC's, we propose that the standard 
scaling behavior of OP can be generalized as O] : 



OP(i, r, L x , L y ) = b~ p / v ^OP**{b~ z t, b^^r, 

b- v w/^L y ,b- l L x ), (8) 

where OP** is a scaling function. Taking b « t x l z in 
Eq.(8) at criticality (t = 0), it follows that: 



OP{t) <xt- p ' Ul - z . 



(9) 



Also taking the derivative of Eq.(8) with respect to r it 
follows that 

dOP(t,T)\ T=0 oct c± ,c ± = — — . (10) 

v±_ z v±z 

Figure 2(a) shows log-log plots of OP(t) vs t, for the IKLS 
model obtained for T = 3.200. Notice that starting from 
a GSC, the power law decay of the OP as described by 
Eq.(9) is obtained after a long time, e.g. for t > 10 5 in 
Figure 2(a). A detailed view of this behavior is shown 
in Figure 2(b) for the four cases studied. Also, fitting 
these curves using Eq.(9) the slopes given by (3/v±z can 
be obtained, as listed in Table I. 

It should be noticed that tiny downward and upward 
deviations from linearity are observed for temperatures 
slightly out of criticality (not shown here). However, us- 
ing GSC, this method for the location of T c is roughly 
twice less sensitive than that used starting with FDC. 
So, our second (independent) estimation of the critical 
temperature is T c /Ki5 = 3.20 ± 0.02. This value is fully 
consistent with our previous estimation made starting 
from FDC, as expected for second-order transitions. The 
same behavior is observed in the other cases studied. 
Data collapsing of the raw data as observed in Figure 
2(a) supports the assumption, implicit in the proposed 
scaling anzats (Eqs. (8-9)), that starting from a GSC 
the only relevant length of the lattice is that perpen- 
dicular to the drive along which a diffusion-like process 
causes relaxation to the NESS. Our results confirm that 
the dynamics does not reveal the anisotropy of the sys- 
tem through the shape of the lattice S — L^^ 11 / L x , in 
contrast to the scalingtreatment of data obtained un- 
der NESS conditions |10|Jll]Jl7|| . Instead, the anisotropy 
enters through the different correlation length exponents 
used in the scaling relationships, i.e., Eqs. (6) and (9). 

Log-log plots of the derivative of Eq.(8) (not shown 
here) support the scaling anzats of Eq.(lO) and allow the 
evaluation of cj_, as listed in Table I. 

Combining the exponents (3/v^z and c± obtained 
starting with GSC it is possible to evaluate 0. Also, 
using those exponents, as well as c 2 and cm as evaluated 
starting from FDC, it is possible to calculate v±, v\\ and z 
(see Table I). Furthermore, for anisotropic systems such 
as the KLS model, the following modified hyperscaling 
relation is expected to hold JlTj : 



+ (d - IK - 2/3 = 7 , 



(11) 



where 7 is the exponent of the susceptibility. So, using 
the numerically evaluated exponents, our estimations of 
7 made with the aid of Eq.(ll) are also listed in Table I. 
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FIG. 2. Log-log plots of the order parameter vs t as ob- 
tained starting with GSC. a) Results corresponding to the 
IKLS model using lattices of different sizes as shown in the 
figure, b) Results obtained for the different models studied 
as identified in the figure. Data obtained using lattices of 
size Lix = 480, L y = 240, except for the FKLS model with 
L x = 960, L y = 60. 



Based on the results summarized in Table I, it is con- 
cluded that: i) All critical exponents can be obtained 
sclf-consistently within the dynamics framework. In con- 
trast to previous numerical estimations performed under 
NESS conditions, the values of the exponents are not- 
biased by any assumption on the shape of the lattice 
used, ii) All the evaluated exponents are consistent with 
Eq.(2), and the predictions of Eq.(l) are far outside our 
error bars, iii) Since two models with macroscopic cur- 
rent (IKLS, FKLS) and another two without it (IRKLS, 
IOKLS), have the same critical exponents, it follows that 
such current is not relevant neither to determine the uni- 
versality class nor for the emergence of the anisotropy. iv) 
The short-time dynamics of the studied non-equilibrium 
systems exhibits universal critical behavior, as already 
observed in equilibrium systems. 
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TABLE I. List of critical temperatures and exponents. The 
symbols (*), (+) and (*+) stand for exponents obtained us- 
ing data corresponding to FDC and GSC, and combinations 
of both, respectively. The predictions of both Eq. (1) and 
Eq.(2), up to the first term of the e expansion, are also listed. 
The results of the oscillatory models are less sensitive to tiny 
changes of T, causing the exponents to be less accurate than 
in the other cases. Notice that z is renormalizated with re- 
spect to the perpendicular direction. 



MODEL 


T c 


<*(*) c,|(*) 0/v±z{+) c±(+) /?(+) 


v±(*+) f||(*+) z(*+) 7 


IKLS 


3.200(10) 0.114(5) 0.406(10) 0.254(15) 0.516(15) 0.330(30) 0.644(40) 1.221(40) 2.016(40) 1.21(11) 


FKLS 


2.780(15) 0.108(5) 0.409(10) 0.245(15) 0.487(15) 0.335(30) 0.668(40) 1.198(40) 2.041(40) 1.20(11) 


IRKLS 


3.160(20) 0.115(5) 0.421(10) 0.228(15) 0.506(15) 0.311(30) 0.671(40) 1.168(40) 2.033(40) 1.22(11) 


IOKLS 


3.160(20) 0.116(5) 0.424(10) 0.235(15) 0.505(15) 0.281(30) 0.635(40) 1.110(40) 2.126(40) 1.18(11) 


Eq. (1) 




wl/8 wl/2 «3/4 w3/4 1/2 


wl/2 «3/2 w4/3 wl 


Eq. (2) 




w 0.114 « 0.410 « 0.262 « 0.532 w 0.33 


w0.63 ss 1.22 « 1.998 « 1.17 
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